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The phonon thermal conductivity of a multilayer is calculated for transport perpendicular to 
the layers. There is a cross over between particle transport for thick layers to wave transport for 
thin layers. The calculations shows that the conductivity has a minimum value for a layer thickness 
somewhat smaller then the mean free path of the phonons. 



The thermal conductivity is a fundamental transport 
parameter Q]. There has been much recent interest in 
the thermal conductivity of semiconductor superlattices 
due to their possible applications in a variety of devices. 
Efficient solid state refrigeration requires a low thermal 
conductivity 0]. Preliminary experimental and theo- 
retical work suggests that the thermal conductivity of 
superlattices is quite low, both for transport along the 
planes HQ, or perpendicular to the planes The 
heat is carried by excitations such as phonons and elec- 
trons. Most theories use a Boltzmann equation which 
treats the excitations as particles and ignores wave inter- 
ference |7]|| . These theories all predict that the thermal 
conductivity perpendicular to the layers decreases as the 
layer spacing is reduced in the superlattice. The correct 
description using the Boltmann equation would be to use 
the phonon states of the superlattice as an input to the 
scattering, but this has not yet been done by anyone. 

We present calculations of the thermal conductivity 
perpendicular to the layers which includes the wave in- 
terference of the superlattice. These calculations, in one, 
two, and three dimensions always predict that the ther- 
mal conductivity increases as the layer spacing is reduced 
in the superlattice. This behavior is shown to be caused 
by band folding in the superlattice. It is a general fea- 
ture which should be true in all cases. The particle and 
wave calculations are in direct disagreement on the be- 
havior of the thermal conductivity with decreasing layer 
spacing. This disagreement is resolved by calculations 
which include the mean- free-path (mfp) of the phonons. 
For layers thinner than the mfp, the wave theory applies. 
For layers thicker than the mfp the particle theory ap- 
plies. The combined theory predicts a minimum in the 
thermal conductivity, as a function of layer spacing. The 
thickness of the layers for minimum thermal conductivity 
depends upon the average mfp, and is therefore temper- 
ature dependent. 

The particle theories use the interface boundary resis- 
tance jl0| as the important feature of a superlattice. A 
superlattice with alternating layers has a thermal resis- 
tance for one repeat unit oiRsL = L\/ K1+L2/ K2+2R B , 
where (Lj,Kj) are the thickness and thermal conductiv- 
ity of the individual layers, and R B is the thermal bound- 
ary resistance. For simplicity assume that L\ = L2 = L, 



which is often the case experimentally. The effective ther- 
mal conductivity of the superlattice is then 
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This classical prediction is that the thermal conductivity 
decreases as the layer thickness L decreases ||. 

The wave theory calculates the actual phonon modes 
uo\{k.) of the superlattice, where A is the band index. 
They are used to calculate the thermal conductivity from 
the usual formula in d-dimensions jl| , 

K <n = E/ ^^ck)K(k)|/ A (k)^^ ( 2) 

where n(u>,T) is the Bosc-Einstein distribution func- 
tion. A rigorous treatment uses Boltzmann theory ap- 
plied to the transport in minibands to find the mean-free- 
path £\(k). At high temperatures, one can approximate 
n ~ k B T '/Tko\(}s) which gives the simpler formula 
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The above formula is quite general. There are two impor- 
tant special cases of constant relaxation time (K T ) and 
constant mfp (Kg) 



(4) 
(5) 



Both of these formulas can be related to the distribution 
P(v z ) of phonon velocities perpendicular to the layers 
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Kg = k B l { dv z P(v z )v z 
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Wave interference leads to band folding jlT 12|. Band 
folding leads to a reduction of the phonon velocities. 
Both K T and Kg are reduced by wave interference. The 
case of constant mfp in one dimension can be reduced to 
a simple formula 
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The integration is eliminated. 

These assertions are best illustrated in one dimen- 
sion, that is with the atomic chain. The simplest model 
for a superlattice has all spring constants identical, and 
the layers differ in their masses. Layer one has N/2 
atoms of mass mi and layer two has N/2 atoms of mass 
777.2 = mi /a. The characteristic matrix for phonons is 
given below for the case N = 6 
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Setting the determinant ||M — w 2 l|| to zero gives the 
characteristic equation (the easiest way to derive it is 
to utilize the similarity of the present problem to the 
Schroedinger equation with the Kronig-Penney potential) 
k(u) of pMll 



cos(kN) = cos(k 1 N/2)cos(k 2 N/2) - 
l-cos(fci)cos(fc 2 ) s . 
sm(fci) sm(fc 2 ) 
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where cos(fci) = 1 — us 2 /2, cos(fc2) = 1 — ui 2 /(2a) de- 
fine the wave vectors {k\,k2) in the individual layers 
in dimensionless units. A typical spectrum is shown 
in Fig.l, where a — 2 and the superlattice periods are 
N = 2, 4, 8, 16. Modes with frequency u> > 2 are confined 
within the layer of the lighter atom, and contribute little 
to the thermal conduction. As the value of N is increased 
in Fig. 1, there is more band folding, and the average 
velocity decreases. The lower curve in Fig. 2 shows the 
thermal conductivity as a function of superlattice period. 
The result for constant mfp is normalized to fc.s-&^2,max, 
where u>2, max = 2 is the maximum phonon frequency 
in the layer 1 of heavy mass. The heat conduction is 
highest at small values of N, and rapidly decreases as 
N increases. It also shows an irregular nonmonotonous 
behavior which we understand and will be discussed el- 
severe. This result is for one dimension. Similar curves 
are found for every case which we have calculated: for 



different values of mass ratio a, and for both Kt and 
K T . Generally, increasing N: (i) increases the amount 
of band folding, (ii) decreases the average velocity in the 
superlattice, and (iii) decreases the thermal conductivity. 

Lower curves in Figs. 3 and 4 show similar calculations 
of Ki for mass ratio a — 2 in two and three dimen- 
sions. Including only nearest-neighbor interactions, the 
characteristic matrix (|l2|) is changed only in its diago- 
nal element, where 2 is replaced by 2 [2 — cos(g)] in two 
dimensions, and 2 [3 — cos(g x ) — cos(q y )] in three dimen- 
sions. Here (q,q x ,q y ) are the wave vectors within the 
planes. These variables are integrated to find the result 
for the heat conduction. These cases also have the fea- 
ture the thermal conducivity falls with increasing value 
of layer thickness. The curve is now smooth, due to the 
averaging over the parallel wave vectors. Similar results 
are found for all values of a. 

In these calculations the thermal boundary resistance 
seems to have disappeared. Note that our model does 
predict thermal boundary resistance for a single inter- 
face [jl6|. The scattering of phonon waves at the inter- 
face causes the thermal resistance. This scattering is in- 
cluded in the present calculation. However, it also causes 
band folding, which is a bigger effect than the thermal 
boundary resistance. The wave calculation predicts that 
in three dimensions the thermal conductivity increases 
as the layer spacing is decreased, which is exactly the 
opposite of the particle calculation which predicts that 
it decreases. Since the effect is due simply to band fold- 
ing, which is a well documented phenomena, then the 
prediction must be accurate. Several prior calculations 
predicted behavior similar to that shown in figures 2,3,4 
|l7],[l8]]. However, no explanation was given for the be- 
havior. 

The missing ingredient in these calculations is the 
mean-free-path of the phonons. When the layer thick- 
ness exceeds the mfp, then interference effects should 
diminish, and the particle model should become appli- 
cable. Our intuition is that the wave model should apply 
when L < £ and the particle model should apply when 
L > i. An phcnomenological method of including I is to 
add a complex part to the wave vector k which is i/£. 
Then recalculate the properties of the superlattice us- 
ing eqn. (|i"3|) . This idea came from Pendry Jljj] who did 
the same thing for electron energy bands (and provides a 
more complete justification for this model): energy gaps 
go away if one includes a finite mfp. A similar result is 
found for phonons. The band gaps diminishc to zero as 
the mfp is decreased for a fixed value of N. 

Upper curves in Figs. 2,3,4 show the thermal conduc- 
tivity in one, two and three dimensions as a function of 
superlattice period for four different values of mfp and 
a = 2. The mfp is given in terms of the number of lat- 
tice spacings (note: not superlattice spacings). Similar 
curves are found for other values of a. For large values of 
I the results are identical to the lower curve (which has 
I = oo). For small values of mfp (£ = 10) the thermal 
conductivity is nearly independent of superlattice period, 
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except where it falls at small values of N. 

The interesting cases have I = 50 or 100. Here the 
thermal conductivity falls as N increases, reaches a min- 
imum, and then starts to increase. This latter behavior 
is the situation expected in the experiments. At room 
temperature, in most solids, anharmonic scattering lim- 
its the phonon mpf to value in the range of 10-100 lat- 
tice constants, which is also the typical value of superlat- 
tice parameter in current devices. Therefore we expect 
the experimental thermal conductivities to behave as the 
curves marked £ — 50, 100 in Fig. 2, 3, 4. The thermal 
conductivity should have a minimum value when plot- 
ted vs. superlattice period. The minimum occurs at 
the cross over between the particle and wave-interference 
types of transport. One experimental result has this be- 
havior |c}|. 
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1. Spectrum of N = 2, 4, 8, 16 superlattices with mass 

ratio a = 2 in the extended zone representation. 
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FIG. 2. Heat conductivity in one dimension as a function 
of superlattice period for mass ratio a = 2 for different values 
of the phonon mfp which is given in units of lattice periods. 
Dimensionless units found by dividing (11) by 2kst- 
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FIG. 3. Thermal conductivity in 2 dimensions as a function 
of superlattice period for mass ratio a — 2 for different values 
of the phonon mfp which is given in units of lattice periods. 
(In units of thermal conductivity of a uniform system of the 
heavier atoms.) 
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FIG. 4. Thermal conductivity in 3 dimensions as a function 
of superlattice period for mass ratio a = 2 for different values 
of the phonon mfp which is given in units of lattice periods. 
(In units of thermal conductivity of a uniform system of the 
heavier atoms.) 
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